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1. Introduction 

The coupling of several angular momenta is one of the important technical problems in 
quantum mechanics. Every standard textbook of quantum mechanics discusses at least 
the coupling of two angular momenta. This problem is of great importance not only 
for practical purposes when dealing with several particles such as in atomic physics, 
nuclear physics, and so on, but is also of quite some mathematical interest. To be 
specific, the problem can be rephrased as a well-known problem in group representation 
theory, namely to obtain a direct sum of irreducible representations of the rotation 
group. There is a large body of literature in which this problem is studied both 
analytically and numerically; see, for example, [HE]- 

The complete analytical formula is known for the coupling of any two angular 
momenta in terms of the Clebsh-Gordan (CG) coefficients. These coefficients 
immediately allow us to express the bases in irreducible representations as super- 
positions of the bases in direct product representations. It is a well established strategy 
to add more than two angular momenta by successive applications of the addition of 
two with the aid of CG coefficients. The addition of many angular momenta can be 
carried out by using the nj symbols, the Racah coefficients, and so on, which are 
generally studied within the recoupling theory. Indeed, the study of many angular 
momenta turns out to be a beautiful mathematical physics problem in its own right. 
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Recently, the study of irreducible representations of many identical spin-j systems, 
or many d-dimensional quantum systems (qudits), has been revived in the context of 
quantum information and quantum computing. When coupling many spin-j systems, 
the structure of the Hilbert space can be described by representation theory of the 
symmetric group. This fact is known as the Weyl-Schur duality which can be utilized 
to solve many problems [3]. Interesting examples of this kind include an estimation 
for a spectrum of an unknown quantum state [4], quantum communication without 
sharing a reference frame [5], universal coding for a classical-quantum channel [6], 
and others [7] . Another interesting feature among these studies is the proposal for an 
efficient quantum circuit to obtain the irreducible representation of the ./V-fold tensor 
product of a d-dimensional Hilbert space, which requires only a total number of gates 
of order TV log(d, log N, log 1/e) up to accuracy e [51 liJl [TU] . 

Let us look at the coupling of several angular momenta using the CG coefficients. 
The first step is to add two angular momenta which are conveniently chosen from all 
angular momenta. The next step is then to add each of the obtained angular momenta 
and another one chosen from the yet-uncoupled angular momenta. By repeating this 
binary coupling many times, one can arrive at the desired result. It is a rather 
straightforward task to perform each step, but the final result cannot be obtained 
in a simple manner. The major obstacle is that the computational complexity of such 
a coupling of several angular momenta grows rather rapidly with the total number of 
angular momenta and the dimension of each angular momentum. 

As the simplest case, we consider the addition of N spin-1/2 systems which can 
be decomposed into a direct sum of irreducible representations labeled with angular 
momentum j as 

^f/2=0^ (1) 
j£J 

where the index set J = {N/2, N/2 - 1, . . .} has (N + 1) /2 or N/2 + 1 elements if N 
is odd or even, respectively. The multiplicity of each irreducible space is 

c . = iV! ^' + 1 ) ^ (2 ) 

3 (N/2+j + l)\(N/2-j)V [) 

When another spin-1/2 system is to be added to the obtained result, we need to couple 
this new spin-1/2 state to J2j c j ~ (2iV)!/(A^!) 2 ~ 4 N different angular momenta, 
which number grows exponentially for large N. 

Another disadvantage of the binary coupling is that the various constituent 
angular momenta are not treated on equal footing. In other words, the resulting 
angular momentum states depend on the way one chooses the paring in the 
intermediate steps. The binary coupling might not be a wise choice when dealing with 
many identical angular momenta. To overcome this problem, a novel coupling scheme 
was proposed more than four decades ago independently by Chakrabarti [11] and 
Levy-Leblond and Levy-Nahas [12] . They studied the non-binary couplings of three 
angular momenta without employing the binary coupling. Their coupling scheme is 
generally referred to as the symmetric coupling or the democratic coupling, which 
reflect the fact that their choice of the complete set of commuting operators (CSCO) 
contains all three angular momenta with equal weights. To our knowledge, there has 
been no generalization of their non-binary coupling to the case of more than three 
angular momenta. 

It is our main motivation here to clarify the meaning of the symmetric coupling 
for four identical spin systems and then to provide a possible solution for the coupling 



Symmetric coupling of four spin- 1/2 systems 



3 



of four spin- 1/2 systems. While the main result was already reported in [13] . we 
have not clarified the meaning of the symmetric coupling as yet. We hope that our 
construction paves the way toward establishing a non-binary coupling of many angular 
momenta. 

This paper is organized as follows. In section [2] we give the mathematical 
background of the Weyl-Schur duality, a possible definition of the symmetric coupling 
of many identical angular momentum systems, and discuss the relation to the 
reference-frame-free subsystems. We then study the symmetric coupling of many 
spin- 1/2 systems in the second largest angular momentum subspace in section [3] The 
detailed analysis in the case of four spin-1/2 systems is shown in section [4] We close 
with a summary and discussion in section [5] 



2. Symmetric coupling of TV spin-1/2 systems 

We briefly summarize some of relevant mathematical facts. Readers are referred to 
[3] for more concise discussions. In the rest of paper, we mainly consider TV spin-1/2 
systems (2-dimensional systems) unless stated explicitly. 



2.1. Weyl-Schur duality 

The TV-fold tensor product of the 2-dimensional Hilbert space can be decomposed into 
the following direct sum structure: 

(c 2 )® N = s v ®n„, (3) 

i/gPar(N,2) 

where Par(N, 2) stands for the partition of TV into two non-negative and non- 
increasing integers, i.e., Par(N, 2) = {{v\, v 2 ) G Z 2 |j/i + v 2 = TV, v\ > v 2 > 0}. In the 
decomposition ©, known as the Wedderburn decomposition, the subspaces 1Z V arc 
the representation spaces for the general matrix group over the complex field with the 
dimension 

r{v) = vi - v 2 + 1, (4) 
and S u are the representation spaces for the symmetric group Sn with the dimension 
TV! { Vl -v 2 + l) 

= —7 — 771 j— • 5) 

[yi + 1)\ v 2 \ 

This dimension s(y) is same as the multiplicity Cj in ([5]) for j = (y\ — v 2 )/2. In other 
words, we can also label the subspaces with a single quantum number j in accordance 
with 

TV . TV . 

v 2 = — -j, (6) 

with j € J. Note that r{v) = 2j + 1 is the dimension of the subspace T>j with angular 
momentum j. As an example, consider the coupling of three spin-1/2 systems, in which 
case the partition is Par(3, 2) = {(3, 0), (2, 1)}. The dimensions of the corresponding 
subspaces are r(3, 0) = 4, s(3, 0) = 1, r(2, 1) = 2, and s(2, 1) = 2. 

Because of the decomposition ([3]), the TV-fold tensor product of the two- 
dimensional non-singular matrices A e GL(2,C), i.e., A® N , acts irreducibly on the 
subspaces 1Z„, and the unitary representations of the permutation operators Pi^...^ 
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act irreducibly on the subspaces S v . Throughout the paper, we denote the permutation 
from 1,2, . . .,N to iii 2 . ..i N by P ili2 ...i N , that is, 

P =( 1 2 N ) (7) 

\ i\ %2 ■ ■ ■ in J 

Therefore, A® N and Pi 1 i 2 ...i JV are decomposed as follows: 

i>£Par(N,2) 

Piii2---iN = S v <E) I r {u)- (8) 

i>£Par(N,2) 

These decompositions are the essence of the Weyl-Schur duality which states that 
operators commuting with all elements of A® N are expressed as linear combinations 
of the unitary permutation operators Pi t i 2 ...i N with complex coefficients. Moreover, 
its inverse also holds, that is, if operators commute with all elements of Pi^.^in , they 
are a linear combination of A® N with complex coefficients. The Weyl-Schur duality 
holds for the general case of TV- fold tensor products of <i-dimensional systems [3]. 



2.2. Complete set of commuting operators and missing label operators 

As a mathematical problem, to calculate the completely reducible representation is 
equivalent to finding the CSCO whose joint eigenstates define the representation 
uniquely up to arbitrary phase factors. A simple counting argument shows that 
the number of elements of the CSCO is N for the case of addition of N arbitrary 
angular momenta^ This follows from the fact that the direct product representation 
is given by the joint eigenvalues of the z-components of the individual angular 
momentum. To be precise, the states are labeled by 2N quantum numbers in 
the direct product representation. Upon denoting the individual spin operator by 
J i = &i/2 {I — 1,2,..., N), where the ats are the Pauli spin operators constituting 
the Lie algebra su(2), the squares are J| = Jg ■ Ji and the z-components of the 
individual angular momenta are 3( z J§| When considering systems in which each angular 
momentum is fixed, we will not write the eigenvalues of the Casimir operator J| 
explicitly. This paper deals with such a system and the total number N of constituents 
thus determines the representation. 

From the general theory of quantum angular momentum, two of the commuting 
operators in the CSCO are immediate. The first one is the Casimir operator of the 
rotational group, which labels the total angular momentum quantum number, and the 
second is the z-component of the total angular momentum. Therefore, the minimal 
number of elements in the CSCO that are still to be constructed is N — 2. These 
remaining N — 2 operators are usually referred to as the missing label operators 
(MLOs), and finding the MLOs has been a standard but rather difficult problem 
in the representation theory [TH [151 US] ■ Note that the CSCO is not unique in general 
and to list all possible families of CSCO seems an untrackable problem except for some 
special cases. 

X The general theorem guarantees that the minimal number of elements of CSCO for finite 
dimensional representations can always be reduced to one. In this paper, however, we adopt the 
direct sum decomposition J3]l and wish to find the CSCO according to this decomposition. 
§ The spin operators act on the full Hilbert space, as illustrated by c?2 = I2 ® 3 §5 I2 ® I2 ® • ■ ■ ® 12- 
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With the total angular momentum J = Y]^—i Je, the above mentioned Casimir 
operator and the z-component of the total angular momentum are J 2 = J ■ J and 
J z = Ji z , respectively. The former operator J 2 specifies the angular momentum 
space j in (fTJ) and the partition v in ([3]). The latter operator determines the 
representation T>j in (fTJ and the representation 1Z V in (|3|). By definition, the MLOs 
commute with J 2 and J z , and their eigenvalues are non-degenerate within each 
subspace specified by the common eigenvalues of J 2 and J z . Without loss of generality, 
these non-degenerate eigenvalues can be chosen as real, and thus the MLOs can be 
given by hermitian operators. It is not difficult to show that the MLOs live only in the 
subspace 5„, and they are expressed as the superpositions of the unitary permutation 
operators. This leads to the key observation that the MLO problem for N spin-j 
systems can be solved by finding a suitable unitary representation of the symmetric 
group S N . 



2.3. Symmetric coupling of three spin- 1/2 systems 

The representation theory of the symmetric group has been much studied [17] . In the 
standard treatment, the irreducible representations can be constructed in real matrix 
forms by employing the canonical subgroup chain 

S N D S N -l => ■■■ => S 2 . (9) 

In the simplest case of S3, for example, there exist three different irreducible 
representations corresponding to three possible Young tableaux. Besides trivial one- 
dimensional representations for the totally symmetric and anti-symmetric subspaces, 
the remaining non-trivial one is the two-dimensional subspace. In the above choice of 
subgroup chain (j9)), one can choose three possible proper subgroups 62(12), 5*2(23), 
and S^^l) where 52(11,12) is the transposition subgroup between two indices (1112), 
e.g., 5 2 (12) = {Pi23,-P2i3}- When the reduction 5 3 D 5 2 (12) is adopted, the 
representations corresponding to the decomposition (|8|) are 

P213 = li U I °, ) i> 



-1 

Pl32 = J l® 7 4©^ (^5 ^f) ®/2, (10) 

and the other representations for the elements of S3 are generated by the combinations 
of these two transpositions. Since a particular proper subgroup 52(12) is diagonalized, 
the symmetry of the subspace representation labeled by the partition v = [y\, v/) is 
determined by the subgroup 52(12), i.e., we have invariance under the transposition 
between the two indices (1,2). 

With these observations, the coupling of three spin-1/2 systems is solved by 
specifying the MLOs. As discussed, the number of MLOs is 1 for the N = 3 
case. This operator is identified with the transposition operator P213 by adopting 
the representation dTUl) . Upon noting that the transposition operator between two 
given systems (k, £) is expressed in terms of the Pauli operators at and ae (k 7^ t) by 

Pkt = iz{h + o t k-b t i )=P lk , (11) 

the MLO is written as F 2 i3 = (^8 + a\ • <?2)/2. It is straightforward to see the 
correspondence with the standard binary coupling, in which the MLO is given by the 
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square of the intermediate coupled angular momenta. In the above choice of MLO 
-P213J the corresponding operator in terms of angular momenta is 

Jh = J12 ■ J12 = 3(3/8 + ffi ■ S 2 ) = P123 + P213, (12) 

where Jjk = Jk + Je is the intermediate angular momentum. Because the identity 
P123 = Is is irrelevant as far as MLOs are concerned, P213 and Jf 2 are essentially the 
same MLO. This illustrates the claim that the coupling of identical angular momenta 
is obtained by the representation theory of the symmetric group. The other possible 
MLOs are also found to be J| 3 = P123 + P132 and J 31 = P123 + ^321- Importantly, 
all three representations are related through the action of unitary transformations, 
and the subject matter of recoupling theory is to study the relationships among these 
different representations. 

We now show a different coupling scheme by reconsidering the subgroup chain 
(0). It is well-known that the cyclic permutation is also a proper subgroup of Sn. 
Denoting the iV-cyclic permutation group by Cat, the alternative subgroup chain is 

S N DC N DC N -!D ...DC 2 = S 2 . (13) 

In the case of three spin- 1/2 constituents, the cyclic group of order 3 is C3 = 
{P123, P231, P312}, and the natural representation of C3 is the diagonal matrix where 
the elements are powers of the basic 3rd root of unity. Thus, two-dimensional 
representations for the elements of S3 can be generated by the combination of the 
3-cycle permutation P231 and 2-cycle permutation P213. They are 

u> 3 

4. 



P231 = h ® h e ( * ^ I ® h 



P213 = h ® h e ( J J J ® h, (14) 

where Ud — exp(27ri/ii) is the basic dth root of unity. One can, of course, check that 
the two representations in (|10|) and (|14j) are related by the unitary matrix 

h, (15) 



T T 1 ( 1 1 

V2 V- 1 1 



but this representation has higher symmetry than the previous one. This is because 
the cyclic permutation subgroup C3 has order 3 rather than 2 for the transposition 
subgroup. We remark that this choice of irreducible representation is only possible 
with complex numbers. 

The MLOs are readily found through the relation P231 = P132P213 and the Pauli 
operator representation of the transition operator pip . 

P231 = j [h + Si ■ S 2 + <? 2 • S 3 + 1S1 ■ (ct 2 x <7 3 )] . (16) 

Since the other element of the cyclic permutation operator P312 is equally good for 
the MLOs, the two cyclic permutations can be combined to give the real one of the 
MLOs as 

K = ^=( p 23i - P312) = -/=^i • {S% x ff a ) 

SO^WJ^W (17) 
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This form of the MLO is more convenient for us for two reasons. Firstly, the 
eigenvalues of K are ±1 and, secondly, K is projected onto the two-dimensional 
subspace v — (2,1). The corresponding angular momentum states in the j = 1/2 
subspace are 

|1/2,1/2;A) = -^(|100K + |010) W3 2A + |001)), 

11/2,-1/2; A) = -i=(|011) W 3 A + |101) W3 2A + |110», (18) 

where A = 1,2 and we denote by |0) and |1) the kets with m = 1/2 and m = —1/2, 
respectively, for the states of the single spin-1/2 constituents. 

Note that this CSCO {J 2 , J z , K} is indeed identical with that of jTTJ Q2], where 
they proved that this set is the CSCO for the general coupling of three angular 
momenta Jt [t = 1,2,3) with the MLO K = J\ ■ (J 2 x J 3 ). Although this novel 
coupling scheme was discovered almost half a century ago, the generalization to the 
case of more than three angular momenta is — to our knowledge — still an open 
problem. The reason seems that there exist infinitely many different MLOs, which are 
equivalent in general. It is then rather hard to conclude that a particular choice of 
MLOs provides the symmetric coupling without a precise definition of the symmetric 
coupling. In the next subsection, we attempt to give a possible definition of the 
symmetric coupling in the case of identical spin systems. 

2.4- A proposal for symmetric coupling 

We define the symmetric coupling of identical angular momenta, not necessarily for 
the case of spin-1/2 systems, as follows: 

If the MLO projected onto a subspace labeled by the partition v can 

be chosen as a linear combination of cyclic permutation operators 

with the possible maximal order in this subspace, then identical (19) 

angular momenta are said to be coupled symmetrically within the 

subspace. 

The word "symmetric coupling" refers to the fact that this coupling respects the 
cyclic permutation symmetry within each subspace. In this definition, the possible 
maximal order of the cyclic permutation operator still needs to be stated explicitly. 
In this paper, we set the order of the cyclic permutation subgroup to be s(y) + 1 
for the partition v. If the order is less than s(y) + 1, we do not have a symmetric 
coupling in the subspace. Note that the cyclic permutation group is abelian and 
hence all linear combinations are equally good as the MLOs as long as they have 
non-degenerate eigenvalues. Another important consequence of the abelian property 
is that the obtained irreducible basis is invariant under the same cyclic permutation. 

We emphasize that this definition is rather limited since it only applies to the case 
of addition of identical angular momenta for a few angular momentum systems. In 
particular, the number s(u) + 1 becomes greater than N for more than four spin-1/2 
systems and we need to refine the meaning of the possible maximal order of the cyclic 
permutation operator properly. According to the definition (|19p . for example, it follows 
that five spin-1/2 systems cannot be coupled symmetrically in the subspace v = (3, 2) 
{j = 1/2) whose degeneracy is 5. Nevertheless, we will show that a symmetric coupling 
of four spin-1/2 systems is possible. The coupling of more spin-1/2 systems and other 
more general cases are left to future studies. 
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2.5. Relation to reference- frame- free subsystems 

Owing to the Weyl-Schur duality, the role of the permutation symmetry is clear 
in the tensor-product Hilbert space. This leads to the idea of constructing a 
reference-frame-free (RFF) subsystem, or a RFF qudit, in the context of quantum 
information theory. A RFF qudit is a d-dimensional subsystem of a composite 
system, which remains invariant under the same unitary transformation on the 
constituents of the composite system. It is called a rotationally invariant qudit when 
considering the invariant subsystems under the collective rotation rather than the 
general unitary transformations. The RFF subsystems have many applications in 
quantum information and quantum computing [5]. The Werner state for two parties 
is the simplest example |18j . and the generalization to more than two parties have 
been studied in this context [HI [20l [21] . 

When dealing with spin-1/2 systems, any unitary transformation is equivalent 
to a rotation (SU(2) = SO(3)). Then, RFF subsystems can be described by a non- 
negative unit-trace density operator in the subsystems, which commutes with the 
AT- fold tensor product of the rotation Uj = exp(m • Jj). Denoting the density operator 
by p (p > 0, tr{p} = 1), the condition for the RFF subsystem reads 



N 

n 



P, e iH - J 



= 0, (20) 



where J is the total angular momentum operator. From the properties of the Weyl- 
Schur duality, it immediately follows that the only possible form of the density operator 
for the RFF subsystem is 

= Ps{v)®Ir{v)i (21) 
ySPar(N,2) 

which is a linear combination of permutation operators -Pj 1 j 2 ...j 2V . Therefore, the 
construction of RFF subsystems is essentially a problem of revealing the algebraic 
relationship between the permutation operators in the various subspaces, that is, to 
analyze the different choices of MLOs. In fact, the following stronger statement holds: 

All possible MLOs for the subspace v are linear combinations of 
a non-degenerate RFF density operator and the projector onto the 
subspace v, (22) 

apf +PI a (u), 
where a and j3 are coefficients. 

In the Lie algebra theory, the algebra formed by RFF states is known as the (universal) 
enveloping algebra of su(2) [22] . 

There are many ways to describe a quantum state defined in the d-dimensional 
Hilbert space C d . The simplest one is to use the d orthonormal kets \k) (k = 1, 2, . . . , d) 
that are orthonormal and complete, 

d 

(k\k')=5 kk ,, ^|fc)(fc|=7 d , (23) 
fc=i 

and thus form a basis in C d . Correspondingly, the d 2 operators 

Qu=\k){i\ (k,£ = l,2,...,d) (24) 
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that satisfy the closure relation QkiQk't' — &lk'Qkt' i form a basis for the d-dimensional 
operator-algebra space, i.e., the d-dimensional matrix ring over C. We call these d 2 
operators the RFF basis operators. The state of a d-dimensional quantum system can 
be represented as a unit-trace d x d matrix that is semi-definite positive. Another 
possible way is to expand the state in terms of the generators for a SU(d) Lie 
group with real coefficients. The standard Gell-Mann matrices together with the 
semi-definite positivity requirement provide a proper d-level quantum state [531 121] ■ 
The third option is to use the unitary Heisenberg-Weyl-Schwinger (HWS) operator 
basis [25 . The complete set of unitary operators is given by U k V l (k, £ = 1,2, ... ,d), 
where the unitary operators U and V have the period d, i.e., U d = V d = 1, and satisfy 
the Weyl commutation relation U k V e = uj^ M V e U k with uid — exp(27ri/ci) as above. 
From the mathematical point of view, all three operator bases and any other ones are 
equivalent in the sense that we can convert one description into the others through a 
bijective mapping, and an advantage over others may show up depending upon the 
problem of interest. In the following, the first method is mainly considered, and the 
basis operators Q v kt are to be constructed. 

The Weyl-Schur duality and the decomposition in flS) already provide the relation 
between any permutation operator and the basis operators Q v u in the subspace labeled 
by a partition v. Firstly, define the RFF basis operators in the subspace v by 

Q v u = \k){t\®I r{v) , ( 25 ) 

where the operator \k)(£\ lives in the subspace S v . Then, they form the RFF basis 
operators satisfying QkiQ'k>i> — S^^'Sik'Q^,, and the permutation operator can be 
written as 

s{u) 

Pi^-.-iN = PklQkl- (26) 

i>ePar(N,2) k,l=l 

Here, the coefficients p ke are elements of the matrix representation of the symmetric 
group within the subspace v. Using the RFF basis operators, the MLO M v can 
be expressed as a linear combination of diagonal elements of them with different 
coefficients, 

»(") 

M v = J2lxQxx- (27) 

A=l 

In the example of three spin- 1/2 systems, the relation (|26j) for the subgroup chain 
(fT3l) reads 



-Pi 2 3 


= Q?i 0) 


+ Qii' 1) + 


^22 ' 


p231 


= Q& 0) 


+ ^ 3 Q 1 2 i' 1) 


+ w|Q 


P213 


= Q?i 0) 


+ Qi 2 2 1) + 


^521 ' 


Pl32 


,0(3,0) 




+ ulQ 


p521 


= Qi 3 i 0) 


+ ^3^12 


+ ^3<9 



(2,1) 
22 



(2,1) 



(P 3 12) t , 



(2,1) 
21 • 



(28) 

Although these equations seem to be over-determined at first sight, i.e., six equations 
for five variables Q^£, only five equations are actually linearly independent. Note that 
in general there are N\ linear equations for the RFF basis operators of which the 
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total number is s (^) 2 - This number is always less than N\ and, therefore, the 

>/6Par(N,2) 

occurrence of linearly dependent equations is generic. This is due to the fact that the 
entire representation space is not exhausted when considering a problem of N spin-1/2 
systems. The complete set of representations can be obtained for the case of N spin-j 
(j = N/2 + 1) systems, where the following relation holds: 

E ^) 2 =iV!. (29) 

j/£Par(N,N) 

By converting equations (|28|) . we obtain the operators Q v u in terms of the 
permutation operators: 

\ Q22 

gl 1 ) I = ±| lui u> 3 I I P 132 |. (30) 
In terms of the Pauli spin operators, they read 




)(3,0) 
'll 



11, 

~h + g (cti • cr 2 + cr 2 • cr 3 + CT 3 • (T\) 



Q^ 1] 1 1 r 1 



Q 



21) j .4 - ^ (^1 ' ^2 +^2 • CT 3 + CT 3 • ffx) 



22 



± TIT 1 ■ (<?2 x * a) ' 



(O^l'^) = g (CTl • CT 2 + W3CT2 • <? 3 + • ai) . (31) 

With these basis operators Qf-f and Q^' 1 ^, we can express any RFF state as a linear 
combination of them, and this completes the analysis of three spin-1/2 systems. 

3. Symmetric coupling for the second-largest angular momentum 

In this section, we investigate a symmetric coupling of N spin-1/2 systems within 
the second-largest angular momentum subspace |13j . Since the second-largest angular 
momentum has N — 1 components and its multiplicity is cjv/2-1 = N — 1, there are 
(N — l) 2 states in the subsystem in total. The basic ingredient is the grading of the 
N spin-1/2 constituents J e = at/2 (£=1,2,..., N), 

N 

£(A) = £)w#jJ. (32) 

The parameter A takes the values A = 1, 2, . . . , N with modulo N and the A = N case 
reduces to the usual total angular momentum operator S(AT) = J = ^/2- They 

satisfy S(A) = E(iV — A) and the commutation relation 

a • S(A), b • S(A')1 = i(S x b) • S(A + A') (33) 
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for all numerical vectors a and b. With this commutation relation, the E(A)s form a 
so-called graded Lie algebra labeled by the integers A (mod N) . 

Denoting the kets for the single spin-1/2 states with m = 1/2 and m = —1/2 by 
|0) and |1) as before, the ket for the state with maximal values of both j and m, i.e., 
ji = mi = N/2, is 

[Ojv) = \0)® N , (34) 

and successive applications of the lowering operator J_ — J x — \J y yield all states for 
the maximal angular momentum space, 



\3um x ) = J .^ +m ^\ , (J-) jl - mi \0 N ). (35) 
V Jr- Ui - m i)' 

The highest states with j 2 = m 2 = N/2 — 1 for the second-largest angular momentum 
space are given by the action of N — 1 lowering operators 

N 

E_(A) = E x (A)-iE„(A) = 5^w^J<_. (36) 

onto the highest state |0jv) in the largest angular momentum space as 

|i 2 ,j 2 ;A) = ^=E_(A)|0 JV ) (A = 1,2,..., AT -1), (37) 
V N 



and successive applications of J_ give all the remaining states |j 2 , m2; A) with m 2 = 
— j2j — J2 + 1, • • • , J 2 - Since E_ and J_ commute with each other, the resulting states 
are 



0/2 + m 2)! 



2/1 ' E_(A)|j 1 ,m 2 + l), (38) 



Y (ji +m 2 + l)(ii + m 2 ) 
which are orthonormal in the j 2 subsystem, 

(j 2 ,m 2 ; A|j 2 ,m 2 ; A') = £u' , 

X] X! b'2,m 2 ;A)0' 2 ,m 2 ;A| = Jj 2 . (39) 

A=l m 2 = -j2 

Here, ij 2 is the projector onto the angular momentum j 2 subsystem. This projector 
is a polynomial of the Casimir operator J 2 = J • J, 



Their orthogonality is more transparent when the states (|38p are written as 

1 N 

|ja,m3;A) = -=5^|€;ja,m 2 )w^, (41) 



where the ./V states |£; J 2 , m 2 ) form a pyramid, 



(£;j 2 ,m 2 \e';j 2 ,m' 2 ) = S m2m > ( <W + — — — ^— ) . (42) 

J2 + W2 + 1/ 
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From the Weyl-Schur duality, the second-largest angular momentum states | j 2 , TTI2] A) 
correspond to the partition v — (N — 1,1) in the decomposition ©, and they are 
(N — l)-dimensional irreducible representations of the symmetric group S^. 

We note that the discrete Fourier transformation that we chose in (|36[) is just 
one of many possibilities for defining the £^(A)s and thus the kets \j 2 , mi; A). More 
generally, any unitary (N— 1) x (N—l) matrix, with Nth-raw matrix elements Uni = 1, 
can serve in S_(A) = YleLi U\e<J_ . For the specific choice of the discrete Fourier 
matrix, the projectors |j2, m 2 ', X)(j 2 , rn 2 ; A| are invariant under the cyclic permutation 
subgroup of order N generated by P23— ivi< 

Following the construction for the = 3 case, it is natural to look for the 
representation of Sn that possesses a cyclic permutation symmetry of order N. In 
fact, without exploring the representation theory, we can immediately construct the 
(N—l) 2 RFF basis operators from the states (|38|) by tracing over the quantum number 
m 2 , 

Q[ N ^= £ \j 2 , m2 ,X)(j 2 , m2 ,X'\ = ( Q { N -^)\ (43) 

m 2 = -j2 

Indeed, we can check the properties 



J,Q 



(N-1,1) _ 



AA' 



= 0, 



n (N-l,l) n (N-l,l) r n (JV-l,l) (A .s 

The explicit construction in (I38|) of the angular momentum states in the subspace 
with j 2 = N/2 — 1, enables us to express the MLO for the second-largest angular 
momentum subspace as a linear combination of diagonal elements of the RFF basis 
operators with different coefficients. It is then straightforward but rather tedious 
to rewrite it in terms of individual Pauli spin operators. In the next subsection, we 
provide an explicit construction of the MLO for the second-largest angular momentum 
subspace in terms of cyclic permutation operators. 



3.1. Missing label operator for the second-largest angular momentum subspace 

In the construction of the second- largest angular momentum states in (|38[) . the 
cyclic permutation of order N is respected according to the subgroup chain of the 
permutation group, i.e., Sn 3 Cn- Utilizing this fact we now construct the MLO for 
the second-largest angular momentum subspace in terms of the Pauli operators. In the 
following, we restrict ourselves to representations of the permutation group elements 
within the second-largest angular momentum subspace whose dimension is N — 1 . 

The cyclic permutation operator C that transforms the index (1,2,. . . ,N — 1, N) 
to (2,3,..., N, 1) is given by 

C = P23— JVi = Pn jv-i-Pjv-1 iV-2 ■ • • P32P21, (45) 
where, as before, Pij denotes the transposition operator between the two indices (i,j), 

P i3 = -(l^ + ffi-ffj). (46) 
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As stated, we choose a diagonal representation of the above cyclic permutation within 
the second-largest angular momentum subspace as 



C = 



32 



to 







N 







\ 



V 



CO 



N-l 
N 



l N-l 



(47) 



where 



-32 



indicates the restriction to the j 2 subspace. The symmetrically 



constructed second- largest angular momentum states are the eigenstates of cyclic 
permutation operator 

C\j 2 ,m 2 ; A) = \j 2 ,m 2 ; X)uj^. (48) 

This is a consequence of the permutation invariance for the largest angular momentum 
states and the commutation relation C£_(A) = u^£_(A)C. 

It is straightforward to construct N orthogonal projectors by the inverse Fourier 
transform of the permutation operators {C, C 2 , . . . , C^ -1 , C N — I 2 n} 7 



1 N 



uj n O . 



(49) 



fc=i 



They have the following representation within the second-largest angular momentum 
subspace 

/Sxi \ 




5A2 







V 



'AT-1, 



(50) 



(51) 



(52) 



S\N-1 J 

and satisfy the orthogonal relation 

P C (\)P C (\') = 5 XX ,P C (\). 

We thus obtain 

P c (A)|j 2 ,m 2 ; A') = \j 2 , vn 2 \ \)6w- 

These relations are obtained directly from the commutation relation P C (A)E_(A') = 
S_(A')_P C (A — A') and the observation that P c (X)\ji,mi) = \ji,mi)5\N. We remark 
that the rank of the projectors P C (A) (A = 1, 2, . . . , TV — 1) is greater than N — 1 
in general. It follows that they project onto not only the second-largest angular 
momentum subspace but other lower angular momentum subspaces as well. 

With the above result, we can express the MLO for the second-largest angular 
momentum subspace as 



JV-l 

i**=E«r i,i): 

A=l 



J P C (A) (53) 
with N — 1 different coefficients <j^ JV_1 ' 1 \ The choice q^ -1,1 ^ = j 2 + 1 — A reads 



Mi. 



N-l 
k=l 



7 fe = ^E(y- A )<- 



(54) 



A=l 
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The coefficients 7*, have the properties 7jv = and "fN-k — ~7fc ; arL d the MLO is 
stated as the summation of C k — C N ~ k with certain coefficients. 



4. Symmetric coupling for four spin-1/2 systems 

In this section the symmetric coupling for four spin-1/2 systems is accomplished 
along the ideas described in section [2j The standard binary coupling of four angular 
momenta is known under names such as the L-S coupling and the j-j coupling [HE]- In 
the first step, two angular momenta J\ and J2 are coupled and similarly the remaining 
two J3 and J4 are coupled. Then, the newly coupled angular momenta J12 and J34 
are coupled in the last stage. With this particular choice of intermediate angular 
momentum states, the CSCO are 

CSC0 12 | 34 = {J 2 , J„ Ji 2 , J 3 2 4 }- (55) 

Since there are three inequivalent ways of paring four angular momenta as the 
intermediate states, there are two more choices for CSCO, namely CSC0 13 | 2 4 and 

CSCOi4|23- 

From the Weyl-Schur duality, the partition of four into two is Par(4, 2) = {(4, 0), 
(3,1), (2,2)}, and the corresponding dimensions of the subspaces for the symmetric 
group are s{v) = 1, 3, 2, respectively. The total number of basis operators to span the 
RFF subsystems is thus l 2 + 3 2 + 2 2 = 14. As noted earlier, this number is smaller 
than the number of elements in S4, namely 14 < 4! = 24. The cyclic permutation 
subgroup of order 4, 

C4 = {-Pi234> -P2341J -P3412J -P4123}; (56) 

is used to diagonalize the subspace (3, 1) (J = 1), and similarly the cyclic permutation 
subgroup of order 3 

C 3 (123) = {^1234,^2314,^3124}, (57) 

is used for the subspace (2,2) (j — 0). We remark that there are two other cyclic 
permutation subgroup of order 4 and three others of order 3. However, all other 
choices lead essentially to the same result. Following the subgroup chain (fl~3|) together 
with the last transposition permutation subgroup 5*2(12), we obtain the following 
unitary representations for the permutation operators: 

' w 4 \ / u, 2 \ 

^2341 =h®I 5 ®\ Qui ®/ 3 © „ 3 (g)Jl, 

ul V W3 / 



2314 



P2134 = h®h®\ 1 + i 1-i \®h@ <8>Ji, (58) 



and all other elements of the permutation group S4 can be generated with these three 
elements. Using the above choice of the representations, the basis operators Q v u can 
be expressed in terms of the permutation operators in a straightforward, if somewhat 
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tedious, manner. For completeness, we list all relations using the notations where 
I = -Pi234> Pij for the transposition permutations, and 





— -P2341, c\ 


— -P3412, cf 


= -P4123 


c 2 


= -P2413, 


cl 


= P4321, cf 


= -P3142 


c 3 


= ^3421, 


eg 


= ^2143: Cl 


= ^4312 


Di 


= Pl342, 


Di 


= fl423 ; 




D 2 


= -P3241, 


Di 


= P42 13, 




D 3 


= -P2431, 


D\ 


= ^4132, 




D A 


= -P2314, 


Di 


= ^3124, 





(59) 

for cyclic permutations as follows. For the subspaces with v = (4,0) and v = (3, 1), 
we have two different representations in terms of odd and even permutation operators 
as 

3 4 



Q33 ' 



Q ( 2 3A) 



'22 



^+E^ 2 +E(^+a 2 ) 
fe=i 1=1 



fc=i 



i^-^^iE^-a 2 



z=i 



^12+^23+^34 + ^41) 

- g(C 2 + Cf + C 3 + Cf) T ^(Ci - C 3 ), 

^z + ^-cl-cf), 
^OPis + ^-Ci-cv 3 ), 



(or) 



t 1 



^ = (oL' ^) 



±£(-l)'(A + D?) + i(C?-CS 



^13 - P24) - -^(Pl2 - P23 + ^34 - P41) 
- ^(^2 + Cf - C 3 - Cf ), 



3 3 

i(l±i)^c 4 3i A + i(l + i)E^ 



;=i 



4(1 + i)(Pi2 - P34) + (1 - i)(P 23 - P41) 



+ 



(l-i)(C 2 -C 3 ) + (l + i)(C 3 -C 3 ) 



(60) 
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For the v = (2, 2) subspace, we have 

"'~ : (D 1+ D 2 2 +D a + Dl) 



12 

UJ 2X 

+ ^-{Dl+D 2 + Dl + D±), (A = 1,2) 
Q ( i2 2) = (Q ( 2i 2) ) f = ^(C 3 + C| + Pia + P 34 ) 

+ ^(C 1 +C 1 3 + P 1 3 + P 24 ) 

+ ^(C 2 + C 2 3 + P 23 + P4l). (61) 

As we see there are many ways of representing the Q% e s by linear combinations 
of permutation operators, and here we choose the following representation to see 
representations for the RFF basis operators in terms of the Pauli spin operators. 
Define the hermitian operators Aj, Kj, and Lj by 

Af = P2134 ± Pi243, = P3214 ± P1432, ^3 = P1324 ± P42317 

K\ — i(Pi342 — pL42 3 ), K2 = i(p3241 — P4213), 
K3 = i(P2431 — P4132), K4 — i(P2314 — P3124), 

L\ — P2143, L2 = P3412, ^3=P4321, (62) 

where the A^s are combinations of elements of the transposition permutation 
subgroup, the Kjs are those of the cyclic permutation subgroup of order 3, and 
the Lj are those of the 2-cycle permutation operators. Their Pauli spin operator 
representations are 

K\ = -\^2 ■ («? 3 X <? 4 ), K 2 = ■ (<?4 X (?i), 

1 1 

K 3 = —^i ■ (<?i x (7 2 ) , K 4 = ■ (<?2 x 

Li = -(he + <?i ■ V2)(h% + ct 3 • <t 4 ), 
L 2 = ^(he + <?i ■ a 3 )(Ii6 + <?2 • <t 4 ), 

L 3 = i(7 16 + 02 ■ CT 3 )(/i 6 + (?4 • <?i), (63) 
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besides the trivial ones for Af. With these notations, we have 



& 0) 



3 

-ll 16 + l(At + A++A+) 
+ ^-(Li + L 2 + L 3 ), 

Q(3.1) 111 1 



'33 



Q22 X) — ^^16 _ — L2 + L 3 ), 



i[(i + i)^r + (i-iM3 



g& 1, = (g& 1, ) t 
= (qL ' ) , 

V 7 ' T (1^1+^2-1^3 -^4)], 

Qif } 1 1 1 



Q 



I 16 (At + At + At) 

(2,2) 4 10 12 



22 

+ ^(ii + L 2 + L 3 ), 

- (Li + W3L2 + w|L a )] . (64) 

This agrees with the result presented in section IV D of [13.. We remark that the 
RFF basis operators for the subspace (2, 2) (j = 0) can be expressed in terms of 
spin-0 singlet states. 

The MLOs are constructed from the diagonal elements with the proper choice of 
coefficients. The most natural form among many different constructions is 

1 x Qf{ r) + x Q 2 | 1] + (-1) x Q 3 3 3 1} = ^{K-y +K 2 + K 3 + K 4 ) = M j=1 , 

1 x Q ( ?{ 2) + (-1) x Q 2 2 ' 2) = -$={K X - K 2 + K 3 - Ki) = M j=0 , (65) 

vl2 



and hence the CSCO, up to trivial multiplicative factors, are obtained as 

CSCO sym = { J 2 , J z ,M j=u M j=0 }. (66) 

Note that this construction of the MLO for the j = 1 subspace agrees with (|51|) of 
section |3H1 
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The corresponding angular momentum states in the j = 1 subspace are 
|1, 1; A) = i(|1000)c4 + |0100)wf + |0010)cjf + |0001)), 

|1,0;A) = -J=[(|1001)-|0110))K + 1) 

+ (|0101> - |1010))(w 2A + 1) 
+ (|0011)-|1100))(wf + 1)], 

|1,-1;A) = -i(|0111)«J + |1011)w2 x + |H01)w2 A + 11110)) (67) 
with A = 1, 2, 3, and those in the j = subspace are 



|0,0;A) = -^=[(|1001) + |0110))^3 A + (|0101) + |1010))^ 3 2A 



1 r LU 1 

+ (10011) + |1100))] (68) 

with A = 1,2. 

From our construction, the symmetry of CSCO sym is clear. To this end, it is more 
convenient to introduce the following five conjugacy classes of the permutation group 
£4 based upon the disjoint cycles of the subgroups by 

CZ(1 4 ) = {P 1234 }, 

C7(21 2 ) = {P2134, -Pl243, -P3214, -Pl432, P.324, -P423l}: 

CI (31) = {P1342, -Pl423, -P3241, -P4213, -^2431, -P4132, -P2314, -P3124}, 

Cl(i) = {i"2341, -Rll23, -P2413: ^3142, -P3421, -P4312}, 

Cl(2 2 ) = {^3412,^4321,^2143}- (69) 

Then, beside a trivial symmetry under the identity -P1234, Mj=i is invariant under the 
cyclic permutation subgroup of order 4 in ([SB")) . 

c 4 M i= iC4 1 = M j= i with c 4 £ C 4 , (70) 

and Mj =0 is symmetric under the 3-cycle conjugacy class and the two 2-cycle conjugacy 
classes, 

czMj^c^ 1 = M J=0 with c 3 € 01(31), 

cjiMpoc^ 1 = M j=0 with c 2 2 € C7(2 2 ). (71) 

We remark that Mj—i is antisymmetric under P32i4,-Pi423 G C7(21 2 ) and 
-P432i,p2i43 € Cl(2 2 ), and Mj =0 is antisymmetric under the conjugacy classes C7(21 2 ) 
and C7(4). These symmetries should be contrasted with the standard binary coupling 
scheme that gives CSCO12134 or the like. For identical spins, the MLOs of the kind 
J 2 2 , J| 4 are then only invariant under two of the classes in C7(21 2 ) and one of the 
classes in C7(2 2 ). This is our main conclusion: The higher symmetry in the MLOs 
distinguishes this coupling scheme from all other coupling schemes. 
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5. Summary and discussion 

In this paper we have given a detailed analysis of the symmetric coupling of four 
spin-1/2 systems by revealing the symmetry of the MLOs. The result is based on our 
proposal for the symmetric coupling of many identical angular momenta in which the 
Weyl-Schur duality plays a central role. The relation between the MLOs and the RFF 
subsystems is also clarified through the analysis. 

An immediate extension is to provide the MLOs for the coupling of four arbitrary 
angular momenta without using the binary coupling. This long-standing open problem 
can be tackled by first solving the case of four identical angular momentum systems 
and then analyzing the eigenvalues of the obtained MLOs in the non-identical case. 

Another interesting problem is to study the possible coupling scheme for N spin- 
1/2 systems in all subspaces along the ideas presented in this paper. 
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